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ABSTRACT: Penguin, or loop, decays of B mesons induce effective flavor-changing neutral 
currents, which are forbidden at tree level in the Standard Model. These decays give special 
insight into the CKM matrix and are sensitive to non-standard model effects. In this review, we 
give a historical and theoretical introduction to penguins and a description of the various types 
of penguin processes: electromagnetic, electroweak, and gluonic. We review the experimental 
searches for penguin decays, including the measurements of the electromagnetic penguins b — > S7 
and B — + K*j and gluonic penguins B — > K-k, B + — ■» uK + and B — + f] K, and their implications 
for the Standard Model and New Physics. We conclude by exploring the future prospects for 
penguin physics. 
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1 INTRODUCTION 



In the Standard Model, flavor-changing neutral currents (FCNC) are forbidden, 
for example, there is no direct coupling between the b quark and the s or d quarks. 
Effective FCNC are induced by one-loop, or "penguin", diagrams, where a quark 
emits and re-absorbs a W thus changing flavor twice, as in the b — ► t — > s transi- 
tion depicted in Fig. [[]. Penguin decays have become increasingly appreciated in 
recent years. These loop diagrams with their interesting combinations of CKM 
matrix elements give insight into the Standard Model. In addition, they are quite 
sensitive to new physics. 




Figure 1: b — > s loop or "penguin" diagram. Of course, in order to conserve 
energy and momentum, an additional particle is understood to be emitted in the 
transition. 

The weak couplings of the quarks are given by the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix of complex amplitudes: 

(Vy4 Vus Vub \ 
Vcd v cs V cb (1) 
V td V u V tb J 

For the Standard Model with three generations, the CKM matrix can be described 
completely by three Euler-type angles, and a complex phase. In the Wolfenstein 
parameterization ||], the CKM matrix is approximated in terms of four real, 
independent, parameters, (A, A, p, rj), which makes clear the hierarchical structure 
of the elements: 

/ 1 - t;A 2 A AX 3 (p-iri) \ 

^ckm ~ -A 1 - iA 2 AX 2 (2) 

\ AX 3 (1- p-irj) -AX 2 1 / 

Since A = sin#c, the well-known Cabibbo angle, is small (A ~ 0.22), this pa- 
rameterization shows that the off-diagonal elements are small and the diagonal 
elements are close to 1. The complex phase, which may be responsible for CP 
violation in the Standard Model, has been assigned to the corner elements in this 
parameterization. 
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1 INTRODUCTION 



1.1 History of Penguins 

The curious name penguin goes back to a game of darts in a Geneva pub in 
the summer of 1977, involving theorists John Ellis, Mary K. Gaillard, Dimitri 
Nanopoulos and Serge Rudaz (all then at CERN) and experimentalist Melissa 
Franklin (then a Stanford student, now a Harvard professor). Somehow the telling 
of a joke about penguins evolved to the resolution that the loser of the dart game 
would use the word penguin in their next paper. It seems that Rudaz spelled 
Franklin at some point, beating Ellis (otherwise we might now have a detector 
named penguin); sure enough the seminal 1977 paper on loop diagrams in B 
decays Q refers to such diagrams as penguins. This paper contains a whimsical 
acknowledgment to Franklin for "useful discussions" Q. 

Prior to 1975, the loop diagram had been neglected. Vainshtein, Zakharov, and 
Shifman [|| discovered the importance of the penguin diagrams, and suggested 
that penguins were responsible for the enhancement of the AI = 1/2 amplitude 
compared to the AI = 3/2 amplitude in weak K — > irir decays. Penguins were 
considered in the B system by Ellis, et al. ||, and determined to be small 
compared to b —* c amplitudes. However, Guberina, Peccei, and Riickl || later 
pointed out that the penguin b — > s decays could have a rate as large as tree-level 
b — ► u decays. From Eqn. ^, we see that penguin decays involve (V^V^I oc A 2 , 
whereas b — > u transitions involve |Ki6pu<f| °c I n 1979, the role of penguins 
in CP violation was pointed out by Bander, Silverman and Soni |?|] who showed 
that interference between penguin diagrams and tree-level diagrams could give 
large CP asymmetries in B decays. In 1982, Eilam j|] added the gluonic penguin 
b — > sg to the inclusive penguin rate. Later large QCD corrections to the radiative 
penguin b — > S7 ||] were calculated. These corrections increased the predicted 
b — ► S7 rate by a factor of ~ 3, enough to make experimentalists sanguine about 
measuring the rate. The inclusive gluonic penguin b — ► sg* was later clarified 
to include the time-like b — > sqq, and space-like bq — ► sq as well as the light-like 
b — > sg |U]], which increased the gluonic penguin rate ]TT| ], 

In 1993, the CLEO Collaboration published the first evidence for electromag- 
netic penguins (Sec. [O]) in the channel B — > K*^. In 1994 they also measured 
the inclusive b — > 57 rate (Sec. |4.2| ) which was in good agreement with the the- 
oretical predictions. In 1997, after many hints in several experiments, CLEO 
found first evidence for gluonic penguin decays (Sec. ||). 



1 . 2 The Importance Of Penguins 

Although s —* u loop diagrams are important in K decays, those decays are 
typically dominated by large non-perturbative effects. A notable exception is 
K+ — * ix^vv (charge conjugate states are implied throughout this review). This 



decay is expected to be dominated by electroweak penguins (Sec. |2.2[) and could 
eventually provide a measurement of \Vtd\- Penguin processes are also possible in 
c and t decays, but these particles have the CKM- favored decays c — > s and t — » b 
accessible to them. Since the b quark has no kinematically-allowed CKM-favored 
decay (Eqn. |||), the relative importance of the penguin decay is greater. The mass 
of the top quark, the main contributor to the loop, is large, and the coupling of the 
b quark to the t quark, \ Vtb\, is very close to unity, both strengthening the effect 
of the penguin. The b — * s (b — > d) penguin transition is sensitive to \Vt s \ (l^tdl) 
which will be extraordinarily difficult to measure in top decay. Information from 
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the penguin decay will complement information on |Vt s | and \Vtd\ from B s —B s 
and B°-B° mixing. 

Since the Standard Model (SM) loops involve the heaviest known particles 
(t,W,Z), rates for these processes are very sensitive to non-SM extensions with 
heavy charged Higgs or supersymmetric particles. Therefore, measurements of 
loop processes constitute the most sensitive low energy probes for such extensions 
to the Standard Model. 



2 THEORY 



In sections 2.1 through 2.5 we give general descriptions of various kinds of pen- 
guins. Details of the effective Hamiltonian theory are given in section 2.6. Sec- 



tions 2/7 and ^8] discuss two important topics in penguin decay: CP Violation 
and New Physics. 



2.1 Electromagnetic Penguins 

In electromagnetic penguin decays such as b — > s*y, a charged particle emits 
an external real photon (Fig. |2|). The hard photon emitted in these decays is 
an excellent experimental signature. The inclusive rate is dominated by short 
distance (perturbative) interactions and can be reliably predicted. QCD cor- 
rections enhance the rate and have been calculated precisely. Assuming uni- 
tarity of the CKM matrix to constrain \Vt s \ the Standard Model predicts |^] 
B(b — ► 57) = (3.5 ±0.3) x 10~ 4 . Unfortunately, uncertainties in the hadronization 
process limit the ability to predict individual exclusive rates from first principles 
of the theory. Phenomenological predictions range from 1% to 40% Jill for the 
ratio R K * = B{B -»■ K*^)/B{b -> s 7 ). 




Figure 2: Feynman diagram for the electromagnetic penguins b — ► S7 and b — > dj. 
The photon can be emitted from the W (shown) or from any of the quarks. 

The electromagnetic penguin decay b — > d 7 is further suppressed by | V^l 2 /! Vt s \ 2 
and gives an alternative to B°-B° mixing for extracting \Vtd\- Experimentally, 
inclusive b — > d'y has large backgrounds from the dominant b — > 57 decays which 
must be rejected using good particle identification or kinematic separation. 



2.2 Electroweak Penguins 

The decay b — * s£ + l~~ can proceed via an electroweak penguin diagram where 
an emitted virtual photon or Z° produces a pair of leptons (Fig. ||a,b). This 
decay can also proceed via a box diagram (Fig. |3|c). The SM prediction for the 
b — > s£ + l~ decay rate is two orders of magnitude smaller than the b — > S7 rate 

pal pen- 
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2 THEORY 



The rate for b — * svv is enhanced relative to b — ► s£ + £~ primarily due to 
summing the three neutrino flavors. These decays are expected to be dominated 
by the weak penguin, since neutrinos do not couple to photons. The predicted 
rate is only a factor of 10 lower than for b — > 57 plf . Unfortunately, the neutrinos 
escape detection, making this decay mode difficult to observe. 



W-^.^ 7 



-k— a*. 




(a) 



( \ W-l iw+ 
s.d b k — ft- — • s.d b ^ — ft- — i s.d 



t 

(b) 



t 



Figure 3: (a) Photon penguin (b) Z° penguin and (c) box diagrams for the 
electroweak decay b — > {s,d)d. + £~ . The diagrams for b — > {s,d)uP are similar, 
except that (a) does not contribute. 



S.5 Vertical Electroweak Penguins 

Another category of penguin is the so-called vertical or annihilation penguin 
where the penguin loop connects the two quarks in the B meson (Fig. These 
rates are expected to be highly suppressed in the Standard Model since they 
involve a b — > d transition and are suppressed by (/b/wib) 2 ~ 2 x 10~ 3 , where 
Jb is the S-meson decay constant which parameterizes the probability that the 
two quarks in the B meson will "find each other", and is the B meson 
mass. The B — > 77 decay is suppressed relative to b — > S7 by an additional 
ctQED- The B — > £ + £~ decays are helicity-suppressed. Because these decays are 
so suppressed in the Standard Model, they provide a good opportunity to look 



for non-SM effects (Sec. 2^ 




B° 



d 



(b) 



/ 70 



£+ 



Figure 4: Vertical or annihilation penguins: (a) B — > 77 and (b) B — > £^ 



2.4 Gluonic Penguins 

An on- or off-shell gluon can also be emitted from the penguin loop. (Fig. 
While the on-shell b — > sg rate had been calculated to be 0(0.1%), the inclusive 



2.5 Other Contributions to Hadronic Final States 
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Figure 5: Feynman diagram for the gluonic penguin b — > sg* . The gluon can be 
emitted from any of the quark lines and can be on-shell or off-shell. 

on- plus off-shell b — > sg* rate includes contributions from b — > sqq and b — > sgg 
which increase the inclusive rate to 0.5-1% [10, 11]. The b — ► dg* penguin rate is 



smaller by |Vtd/Vt s | . Unfortunately, there are several difficulties associated with 
gluonic penguins. There is no good signature for the inclusive b — > sg* decay, 
unlike the b — > 57 case. The branching fraction of individual exclusive gluonic 
penguin channels is typically quite small and hadronization effects are difficult 
to calculate. In addition, many gluonic penguin final states are accessible via 



other diagrams (Sec. 2.5 ) so the gluonic penguin is difficult to assess. Thus the 
"smoking-gun" penguin processes such as B° — ► (fiK° that have contributions 
only from gluonic penguins are eagerly sought. 

2.5 Other Contributions to Hadronic Final States 

While the gluonic penguin gives rise only to hadronic final states, several other 
processes can contribute to the same final states. One important contribution is 
from the tree-level b — > u decay. For example, the b — > uus transition (Fig. |6|a) 
and the b — ► sg* penguin transition both contribute to B° — > K + ir~. However, 
the b — ► uus transition is Cabibbo-suppressed, so the penguin process is expected 
to dominate. On the other hand, in B — > ir + 7r~ for example, the small b — ► dg* 
contribution is expected to be dominated by the non-Cabibbo-suppressed tree- 
level b — ► uud transition. In general, most decays to hadronic final states with 
<fi mesons or non-zero net strangeness are expected to be dominated by gluonic 
penguins and hadronic final states with zero net strangeness are expected to be 
dominated by tree-level b —* u. 

Electroweak penguins also contribute to hadronic final states. Every gluonic 
penguin diagram can be converted to an electroweak penguin by replacing the 
gluon with a Z° or 7 (Fig. ||b). Electroweak penguins with internal Z° or 7 
emission are suppressed relative to the corresponding strong gluonic penguin. In 
the hairpin process (Fig. ^c) the gluon, Z°, or 7 is emitted externally and subse- 
quently forms a meson (similar to the leptonic electroweak penguins in Figs. ||a 
and b). External gluon emission is OZI-suppressed p2[ : the color-octet gluon 
has difficulties forming a color-singlet meson! These hairpin processes, such as 
B — ► (fin, are expected to be dominated by electroweak penguins. A possible 
exception involves decays such as B — > n'K, where it has been suggested that a 
gluonic-hairpin diagram could be significant (Sec. 5.1.3| ). 



The vertical electroweak penguin diagram, similar to Fig. ||b, with the lepton 
pair replaced by a di-quark pair, is highly suppressed and is important only for 
decays such as B° — ► cficfi, where no other diagrams contribute. 

In the annihilation diagram the b and u quarks in a B~ meson annihilate to 
form a virtual W~ . The annihilation diagram is suppressed by \V u f,\ and by 
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2 THEORY 



JbI^b and is expected to be mostly negligible. In the exchange diagram, a 
b —* u transition and a d — > u transition occur simultaneously via the exchange 
of a between the b and d quarks in a B° meson. The exchange process is also 
suppressed by \V u b\ and fs/TriB, and is also expected to be negligible, except in 
decays such as B° — ► K + K~ where no favored diagrams contribute. 



b — »" — *~ " » w 

(a) (b) (c) 

Figure 6: Examples of other diagrams which contribute to hadronic final states: 
(a) tree- level Cabibbo-suppressed b — ► ims, (b) electroweak penguin (c) hairpin 
diagram. 




2. 6 Effective Theory of Penguin Decays 

At high energy scales, \x ~ M\y ~ 80 GeV, quark decays are governed by Feynman 
diagrams such as those depicted in the previous sections. To obtain an effective 
low energy theory relevant for scales /i ~ m& ~ 5 GeV, heavy degrees of freedom 
must be integrated out to obtain an effective coupling for point-like interactions 
of initial and final state particles p^] . For semileptonic decays (e.g. the familiar 
(3— decays in nuclear physics), this integration corresponds to derivation of the 
Fermi theory of point-like four-fermion interactions from Electroweak Quantum 
Field Theory. The effective theory relevant for penguin decays is obtained by 
generalization of the Fermi theory, as depicted in Fig. 0. The heavy degrees of 
freedom in loop decays are W, Z° and t. After the integration they don't appear 
explicitly in the theory, but their effects are hidden in the effective gauge coupling 
constants, running masses and, most importantly, in the so-called Wilson coeffi- 
cients (Cj) describing the effective strength of the local operators (Qi) generated 
by electroweak and strong interactions. The operators can be grouped into three 
categories: i = 1, 2 — current-current operators (Fig. ||a); i = 3, . . . , 6 — gluonic 
penguin operators (related to diagrams in Fig. ||); i = 7, . . . , 10 — electro- weak 
penguin operators (Figs. || and ||). The effective Hamiltonian for b — ► s penguin 
decays has the following form: 

AC 10 

n eS = — /=f v tb v t * s Cibi) Q l (u) 

V 2 i=l 

Technically, the calculations are performed at a high energy scale u ~ Mw, and 
then evolved to a low energy scale u ~ nif, using renormalization group equa- 
tions. This evolution mixes the operators: Ci(u) = ^ Uij(u, Mw)Cj(Mw). The 
renormalization guarantees that the u dependence of Cj is canceled by the /i 
dependence of Qi, thus any observable quantity should not depend on the renor- 
malization scale u. However, since the calculations are performed perturbatively, 



2.6 Effective Theory of Penguin Decays 
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the truncation of the perturbative series induces [i dependence of the theoreti- 
cal predictions, which often dominates the theoretical uncertainty. Higher order 
terms must be included to minimize the \x dependence. 







1,9 ^1,9 



Figure 7: Derivation of effective low energy theory from high energy Quantum 
Field Theory. Derivation of the Fermi theory of the W exchange diagram for 
semileptonic b — > ulv decay decay") is shown at the top. Generalization for 



the loop processes b — ► 57 and b — > sg are shown at the bottom [24]. 



Even though renormalization mixes the operators, specific processes are mostly 
sensitive to a small subset of Wilson coefficients. For example: b — > 57 to C7, 
b — ► sl + l~ to C7, C9 and C\q. 

Since extensions of the Standard Model contribute additional diagrams at the 
high energy scale, they modify the values of the Wilson coefficients in the effective 
low energy theory. 

While Wilson coefficients represent short distance (i.e. high energy) electro- 
weak and strong interactions, the operator elements < X\Qi\B > are influenced 
by long distance (i.e. low energy) strong interactions (here \B> represents the B 
meson and \X > its decay mode). Therefore, unlike the Wilson coefficients, the 
operator elements cannot be obtained perturbatively due to the confining nature 
of strong interactions at large distances. 

Fortunately, when \X > represents an inclusive final state, expansion in powers 
of 1/mb shows that to leading order < X\Qi\B >«< s\Qi\b>, where <s\Qi\b> is 
an operator element for free quarks which can be easily calculated. The first non- 
perturbative corrections are of second order [25|, 0(l/mft 2 ), and are small, thanks 



to the heavy b quark mass. Experimentalists must sum over all possible hadronic 
final states to determine an inclusive rate. For example, B(B — > X s 7) B(b — ► 
57), where X s represents a collection of charmless hadrons with net strangeness 
-f. 

When \X> is an exclusive final state it is difficult to obtain <X\Qi\B > from 
first principles. Numerical treatment of strong quantum fields (Lattice QCD) has 
been useful so far only for the simplest cases in which part of the final state is 
non-hadronic, \X >= \hL >, such as B — > K*j. Since leptons and photons are 
not involved in long distance interactions, the operator element factorizes into 
hadronic and non-hadronic currents, < hL\Qi\B >=< h\J\\B >< L\ J2IO >. The 
latter can be written explicitly. Non-leptonic final states are the most difficult 
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to calculate. Heavy Quark Effective Theory, which is so useful for describing 
ordinary b — » c decays, is of little use for b — ► s, d decays, since the final state 
quarks are light. Phenomenological models used to predict rates for gluonic 
penguin decays make many assumptions, the accuracy of which is often difficult to 
assess, and the scope of the predictions is usually limited to two-body final states. 
The approach which has often been employed p6|, 27, 28 1 is based on factorization 



[pS|]. The rationale of factorization in hadronic decays lies in the phenomenon 
of color-transparency [|30]|, in which one expects intuitively that a pair of fast- 
moving quarks (in two-body decays Eh ~ mg/2) in a color-singlet state effectively 
decouples from soft gluons. Therefore, long distance final state interactions (FSI) 
can be neglected. Short distance FSI mediated by hard gluon exchanges can be 
included perturbatively. The latter are important for predictions of the strong 
phases which make direct CP violation possible. With the factorization ansatz, 
the matrix elements < hi /12 1 Qi \ B > can be expressed as a product of two hadronic 
currents: <h\\ J\\Bxh2\ J2|0>. This involves both the matrix elements of the 
singlet-singlet and octet-octet currents Jj . Since the octet-octet matrix elements 
are not directly measured they are usually discarded. To compensate for this, the 
effective strengths of the singlet-singlet current matrix elements are renormalized 
by replacing the inverse of the number of colors (1/N C ) by a phenomenological 
color parameter £. While the assumption of factorization works well for tree- level 
b — > c decays g3], its applicability to penguin and b — > u decays needs to be 
verified. Also it is not clear if a single parameter £ will suffice to incorporate 
non-perturbative effects in all ten operators of the effective Hamiltonian. 



2.7 CP Violation 

Unitarity of the CKM matrix (Eqn. |l]) leads to several constraints, the most 
interesting of which is the orthogonality of the first and third columns: 

V u dV* b + V cd V* b + V td V t t = (3) 

Eqn. H defines a triangle in the complex plane (Fig. |8|) called the Unitarity Trian- 
gle. The lengths of the sides of the Unitarity Triangle are given by the magnitudes 
of the CKM matrix elements. The angles are given by the phase of the CKM 
matrix elements: j3 is the phase of V ub , 7 is the phase of V t d, and a = ir — f3 — 7. 
Physicists wish to measure the sides of the triangle, and independently measure 
the angles in order to check the consistency of the Standard Model. 

The CPT theorem forbids partial rate asymmetries which are CP violating to 
occur at tree level. The addition of the penguin diagram, with different CKM 
matrix elements and thus different weak phases, allows quantum interference 
between two amplitudes, and therefore can produce direct CP violation. CP 



2.8 Non-Standard Model Possibilities 



11 



violation caused by interference between the penguin and tree processes can lead, 
for instance, to rate asymmetries in B + vs. B~ decays. There are also schemes 
to measure 7 or a using relations between tree-level and penguin amplitudes and 



isospin relations in B — > Kir and related decays [33|. Many of the most promising 
measurements of CP violation in the B system rely on interference caused by B — 
B° mixing (indirect CP violation). For example, the penguin decays B° — > (fiKg 
and B° — ► rj Kg can be used to measure (3 [|35], |3(|. The tree-dominated decay 
B° — > 7r + 7r _ can be used to measure a; however, additional penguin diagrams can 
cause direct CP violation as well as indirect, which may complicate extraction 
of a. This effect is known as penguin pollution. 

2.8 Non-Standard Model Possibilities 

Penguin diagrams are a very good place to look for new physics at low energies 



[37]. Loops are very sensitive to additional heavy particles, for example, charged 
Higgs, SUSY particles, or fourth generation particles. The inclusive b — ► 57 rate 
measured at CLEO has already placed constraints on charged Higgs models, and 
an anomalous WWj coupling, and other models (Sec. |4.2| ). 

Non-standard models with four generations of quarks can, while respecting 
the experimental measurement of b — ► S7, significantly decrease the b — > d'j 
electromagnetic penguin and b — > d gluonic penguin rates while noticeably in- 
creasing the b — > d electroweak and vertical penguin rates. Models with Z°- 
mediated flavor-changing neutral currents can enhance B° — > £ + £~ by two orders 
of magnitude over the Standard Model value and B® — > £ + £~ by one order of 
magnitude. These models can also increase the electroweak penguin-dominated 
decays B + (B®) — ► </>7r + ((/>7r°) by two (one) orders of magnitude above the SM 
prediction without violating the current limits on b — > s£ + £~ or B° — > Kg Kg. 
Multi-Higgs doublet models can also enhance the rates of di-lepton processes 
such as B° — > £ + £~, and in particular, b — > st + t~ can be enhanced without 
affecting b — ► s/i + /i~ or 6 — > se + e~. In supersymmetric models, contributions 
from a charged Higgs in the penguin loop can be cancelled by contributions from 
charginos and gluinos in the loop, thus leaving b — > S7 at the Standard Model 
value, while increasing the b — > s£ + £~ rate. Similarly, if gluinos couple more 
readily to gluons than to photons, the b — ► sg* rate can be strongly enhanced 
without affecting b — > S7 [38|| . In particularly favorable scenarios, SUSY "pen- 
guino" effects can dominate Standard Model penguin effects in decays such as 
B° — ► K <fr [^9[. Even if new physics conspires to give the same rates as Standard 
Model predictions, there will likely be effects in CP violation. For instance Lon- 
don and Soni |36| point out that differences between the value of sin 2/3 measured 
from B — > tpKg and the value measured from penguins such as B° — ► 77' if§ might 
indicate the presence of non-SM processes. These are just a few of the many ways 
that new physics can be detected in B penguin decays. 

Recently, there has been much discussion about new physics enhancing the b — * 



sg* rate |38|, 4C, 41, |42j], thereby solving a couple of "mysteries" in B physics: the 
B semileptonic branching fraction is measured j£| to be smaller than theoretical 
prediction s jB4j ; and the number of charm particles per B decay is smaller than 
predicted |45|| . Both of these mysteries can be solved by increasing the non- 
leptonic, non-charm B width, e.g., b — * sg* which is so far not experimentally 
well-constrained. A rate of 0(10%) seems sufficient. However, limits on the 
b — ► sg* rate are starting to rule out this explanation (Sec |5~4| ), 
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3 EXPERIMENTAL OVERVIEW 

At present b quark decays are under investigation with high statistics data sam- 
ples from three different colliders: the Cornell Electron-positron Storage Ring 
(CESR) producing BB pairs in decays of the T(45) resonance just above the 
e + e~ — > BB threshold; the Large Electron-Positron collider (LEP) at the Eu- 
ropean Laboratory for Particle Physics (CERN), which produced bb pairs in 
decays; and the Tevatron at Fermi National Accelerator Laboratory (FNAL) pro- 
ducing bb pairs in pp collisions. Various production aspects of these machines are 
compared in Table [fl. 

Table 1: Various parameters characterizing experimental environments at three 
colliders used in analyses of b quark data. Explanation of the symbols: 
a— cross section, C— luminosity, (3— velocity of b quarks, fi^jcr— mean decay path, 
/—fractions of b hadron species produced. 



Quantity 



T(45) (CESR) 



Z° (LEP) 



Tevatron 



cr(66) (nb) 
a(bb)/a(qq) 

C peak ( 1Q 31 clR -2 g -l ) 

J Cdt (pb _1 ) analyzed 
in pipeline 
66 pairs analyzed (10 6 ) 

P 

f3 r jcT (fim) 
fragmentation 

background 
B energy 
spatial separation 

of 6 and 6 
fs+ ~ Jb" 

/b° 

/a 6 



1.1 

- 0.3 
48.0 
3100 
3600 
3.3 
> 0.07 
30 



Ehc 



no 
0.5 



9.2 
> 0.2 
1.1 
120 



0.9 
- 1 
2600 

some 

OJE'beam 

yes 

- 0.4 

- 0.1 

- 0.1 



30000 (~ 6000) t 
- 0.001 
2.5 
- 110 

- 3300 (~ 660)+ 
- 1 
500 

large 



yes 

■ 0.4 

■ 0.1 
- 0.1 



main advantage 



simple production, 
statistics 



vertexing, 
one 6 at a time 



cross-section, 
vertexing 



f These numbers correspond to the central region, \y\ < 1, with high transverse mo- 
mentum of the 6 quark, p t > 6 GeV (the number in parentheses). 



CESR provides the highest luminosity and a simple production mechanism. 
Since no fragmentation particles are produced, the energy of reconstructed B's 
can be constrained to the beam energy which provides for a powerful reconstruc- 
tion technique. Also, decay products from the two B J s in the event populate 
the entire solid angle, allowing background discrimination based on event shape. 
However, since the two B mesons are produced almost at rest, detection of a 
detached secondary B decay vertex is not possible on an event-by-event basis. 

At LEP, since the momentum of the b hadrons is appreciable, their decay 
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products are easily separated into back-to-back hemispheres. A large decay length 
gives rise to many important analysis techniques. Even though the bb cross section 
is much larger than at the T(4S*) resonance, the high energy of the beam limits 
achievable luminosity. The small sample size is the limiting factor for the LEP 
experiments. 

The main asset of the Tevatron experiments is their huge production cross- 
section. A large background cross-section is the main obstacle to overcome in 
these experiments. Specialized triggering is needed to limit data acquisition to 
manageable rates. So far, the results from the Tevatron have been mostly limited 
to channels containing high pt muons. Again, vertexing is an important selection 
tool. 

CESR houses only one experimental apparatus (CLEO II). There are four ex- 
periments at LEP (ALEPH, DELPHI, OPAL, and L3), and two at the Tevatron 
(CDF and DO). All these collider detectors have a similar "onion" structure. The 
beam collision point is surrounded by a thin vacuum pipe, followed by subsequent 
layers of nearly cylindrical detectors. Most of the experiments also have end-caps 
to maximize solid angle coverage. The innermost layer is created by a silicon 
strip vertex detector used to pinpoint production points of charged particles. 
The vertex detector is followed by a larger gaseous charged-particle tracker. To- 
gether they are used to determine particle momenta from curvature in a solenoidal 
magnetic field (except for DO which has no magnetic field). In addition, most 
experiments measure specific ionization {dE/dx) in the tracking devices to obtain 
partial charged hadron identification. To aid particle identification and trigger- 
ing, many detectors also have Time-of-Flight (ToF) scintillation counters sur- 
rounding the tracking system. The DELPHI experiment has a Ring Imaging 
Cerenkov instead, providing superior charged hadron identification. In the next 
layer, an electromagnetic calorimeter measures electron and photon energies by 
integrating over electromagnetic showers developing in a dense medium. Espe- 
cially noteworthy are the CLEO II and L3 scintillating-crystal calorimeters, which 
have superb energy resolution. The final layer comprises a hadronic calorimeter 
which measures energies of charged and neutral hadrons. The calorimeters also 
identify muons which penetrate to the outermost layers. CLEO II has no hadron 
calorimeter; instead a thick layer of iron is used to identify muons. 

In our article, we also refer to the previous generation of experiments at Y(4f>), 
CLEO I at CESR and ARGUS at DORIS (a CESR-like e+e" storage ring which 
operated at DESY in Hamburg), and in pp collisions, UA1 at Super Proton- 
antiproton Synchrotron (SppS) at CERN. The amount of data collected by these 
older experiments was about two orders of magnitude smaller than in the con- 
temporary experiments. 

Future B experiments are discussed in section [6|. 

4 ELECTROMAGNETIC PENGUINS 

4-1 B — > X'-f Exclusive decay modes 

After a b quark decays to S7 via the penguin diagram (Fig. |||) the produced 
s quark and the spectator q (q = u for B~ , and d for B°) turn into hadrons. 
The final state usually contains one kaon and at least one pion (B — > is 
forbidden by angular momentum conservation). Hadronization may proceed via 
creation of an intermediate strange resonance: If* (892), i£"i(1270), -fTi(1400), 



14 



4 ELECTROMAGNETIC PENGUINS 



etc. The existence of penguin decays was first confirmed experimentally by the 
CLEO observation |§ of the exclusive decay B -> K* (892)7, with K* -> Kit. 
The initial observation was based on 1.5 million e + e~ — > T(4S) — > BB events. 
Reconstruction of exclusive final states from B mesons produced at the T{AS) 
benefits from the beam energy constraint: Eb = Et, eam . Thus, energies of the B 
decay products must add up to the beam energy: AE = (Ek* + — Eb eam 0. 
Also the B meson mass resolution is improved by an order of magnitude with 



the use of the beam constraint: Mb = \J Eb eam 2 — (pk* + P-y) 2 - These tight 
kinematic constraints are crucial in background suppression and signal extraction. 
Eleven signal events were observed over a background of two events, estimated 
from the AE and Mb sidebands. 

Since the first observation, CLEO has presented an updated analysis based on 
larger statistics (2.6 million BB events in 2.4 /6 _1 of integrated luminosity) and 
improved analysis techniques [47]. Instead of cutting on various variables to de- 
fine the signal region, and then using sidebands to estimate the background, the 
improved analysis used a maximum likelihood fit to determine signal and back- 
ground yields. In addition to Mb and AE, event shape variables and Mk* also 
were used in the fit. The event shape information was optimized to distinguish 
between the signal and the dominant background due to continuum production 
of lighter quarks (e + e~ — ► qq, q = d,u,s,c). This method improved the signal 
efficiency by a factor of two. Even though the background also increased, the 
signal sensitivity increased by about 30% beyond the gain from the increased in- 
tegrated luminosity. Averaging over various charge modes (K + ir~ , Kg-K , K s tt~ , 
K--k°) CLEO obtained: B(B -► if* 7) = (4.2 ± 0.8 ± 0.6) x 10" 5 . Projections of 
the maximum likelihood fit onto Mb, AE, and Mk* are shown in Fig. |9] for the 
K + tt~ channel which has the largest statistics. 





Figure 9: Projections of the maximum likelihood fit (solid lines) onto Mb, AE, 
and Mk* for the B° — > (K + ir~)'y data (histograms). The horizontal scales are 
in GeV. The vertical scale gives number of events per bin. 



The measured branching fraction is in the range predicted by the Standard 
Model, (1 — 15) x 10~ 5 [16|. It is also an order of magnitude larger than would 



be expected if the penguin diagram were not present [48]. 

LEP experiments looked for these decays in e + e~ — > Z° — > bb data but were 
not able to observe the signal due to an insufficient number of bb pairs. They also 
looked for B s — > ^7 decays. ALEPH set a 90% C.L. upper limit of 29 x 10~ 5 on 
the branching fraction for these decays p9^|. 



4.2 Inclusive measurements of b 



S7 rate 
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4-2 Inclusive measurements ofb — > 57 rate 
4.2.1 Experimental results 

The measured rate for the exclusive mode B — ► if* 7 is in the ball-park of 
Standard Model predictions. Quantitative tests of the Standard Model with 
rates measured for exclusive channels are severely handicapped by our inability 
to calculate hadronization probabilities from the first principles of the theory. 
Predictions of phenomenological models for the K* fraction in b — > s 7 decays, 
R K * = B{B -> K*i)jB{b 87), vary from 1-40% |lfj. One should notice 
improvement in recent lattice-QCD calculations in this area |1S|| . 

Fortunately, when summed over all possible final states, hadronization prob- 
abilities aren't relevant and the inclusively measured rate should reflect the 
short distance interactions which can be accurately predicted using the effective 
Hamiltonian of the Standard Model. Since the first non-perturbative correc- 
tion is expected to be of second order in the Aqcd/ 77 ^ expansion, it should 
be small, thanks to the heavy b quark mass. Next-to-leading order pertur- 
bative calculations have been recently completed for b — > 57. Assuming uni- 
tarity of the CKM matrix to constrain \Vt s \ the Standard Model predicts [O]: 
B(b -> 57) = (3.5 ± 0.3) x 10" 4 . 

When reconstructing simple exclusive final states such as B — ► K* 7, K* — ► 
Kit, backgrounds are usually low due to tight kinematic constraints, in this case 
constraints to the B and K* masses and the beam energy. Inclusive measurements 
are more challenging and they are often background limited. 

The main background for CLEO is again from e + e~ — ► qq events. This back- 
ground can be subtracted reliably with data taken below the e + e~ — > BB thresh- 
old. However, statistical fluctuations in the background level can easily swamp 
the signal unless the background is efficiently suppressed. Backgrounds from B 
decays are less serious since b — > 57 decays are quasi-two-body and produce 
higher energy photons [E~ ~ m&/2) than photons from typical B decay modes. 

CLEO used two complementary approaches to suppress the continuum back- 



ground [50|. In one approach only the photon among b —* S7 decay products 
was explicitly reconstructed. Shape differences between BB events and qq back- 
ground were used for background suppression. BB events are nearly spherical 
since the B mesons are nearly at rest at the T(4S'), while qq events have a distinct 
two-jet appearance. For the best sensitivity all shape variables were combined 
withe the use of a neural-net technique. The signal amplitude was extracted from 
a one-parameter fit to the neural net output variable, with the signal shape and 
the BB backgrounds taken from Monte Carlo simulation, and the continuum 
background subtracted using the below-threshold data. In the second approach, 
all products of the b — > sj decay were reconstructed as in exclusive reconstruc- 
tion. Thus, the constraints to the B mass and beam energy could be used. The 
final state recoiling against the photon, denoted X s , was required to contain a 
kaon candidate (a charged track consistent with by dE/dx and ToF, or a 
Kg — > 7r + 7r~ candidate) and 1-4 pions (including at most one 7T°). In Fig. 10, we 
show the photon energy spectra measured by CLEO with these two methods in 
a sample of 2.2 x 10 6 BB events. The first method has rather large continuum 
background but also high signal efficiency (32%). The second method is very 
good in suppressing continuum background, but the signal efficiency is much 
smaller (9%). The sensitivity of these two approaches is nearly equal, and the 



16 



4 ELECTROMAGNETIC PENGUINS 



measurements of signal amplitudes are only slightly correlated. After combining 
these two methods, CLEO measured B(b -> 57) = (2.32 ± 0.57 ± 0.35) x 10~ 4 , 
in agreement with the Standard Model calculations. 

The X s mass distribution from the CLEO inclusive ^-reconstruction analysis 
is shown in Fig. 12. A clear K* (892) peak is observed followed by a dip which is 
expected since the next excited kaon resonance is l£i(1270). A broad enhance- 
ment at and above the Xi(1270) is observed; this is also expected since many 
resonant states exist in this region. The present experimental statistics are in- 
sufficient to establish a positive signal for any of the resonances beyond K*(892) 
taken separately. 

Combining the inclusive and the exclusive measurements, CLEO determines 
Rk* = (18.1 ± 6.8)% in agreement with some phenomenological estimates [16]. 
In particular the calculations which take the most from QCD, QCD sum rules 
and recent lattice-QCD calculations [19], are in good agreement with the 
data. 

Inclusive b — > 57 decays also have been observed recently in Z° — > bb data by 
ALEPH |j52~|| . The ALEPH analysis, based on 0.8 million bb pairs, searches for 
these decays by selecting events with a single high energy photon and a jet in the 
opposite hemisphere which is 6— like from its observed detached vertex. The B 
hadron candidate in the photon hemisphere was reconstructed by combining the 
photon with other particles in the same hemisphere until the ensemble matches 
the B meson mass within the experimental resolution. The particles were selected 
by their probability of belonging to a B meson which was assigned according to 
their momenta and impact parameters at the interaction point. Up to eight parti- 
cles were allowed in addition to the photon, including charged tracks, 7r° mesons, 
and mesons detected in the hadron calorimeter. Once the B candidate was 
reconstructed, the photon was boosted to its rest frame where the energy E* is 
nearly monochromatic for signal photons (E* ~ rrib/2). Further background sup- 
pression was achieved by imposing requirements on the boosted sphericity and 
photon angle in the rest frame of the B meson candidate. 

To measure the background the sample was divided into eight subsamples, 
seven of which contained little b — > S7 signal. The variables used for this division 
were: lateral shower size in the electromagnetic calorimeter to distinguish prompt 
photons from merged ir° mesons; vertex detachment in the opposite hemisphere 
to distinguish bb events from qq backgrounds; and energy of the B candidate 
to discriminate against final state radiation background. The shapes of the con- 
tributing backgrounds (split into four major components) and of the b — > S7 signal 
are computed from Monte Carlo while their normalizations are fixed by a fit to the 
E* distributions in the eight subsamples. Fig. 11a shows the E* distribution in 
the signal sensitive subsample for both the data and the background and Fig. lib 
shows the data after background subtraction. The total reconstruction efficiency 
(~ 13%) is similar to the efficiency obtained in the CLEO inclusive B reconstruc- 
tion analysis. In spite of 2.75 times fewer bb events, ALEPH is able to observe 
a significant inclusive signal. This should be attributed to a better suppression 
of the light quark backgrounds (the dominant background at the T(4S)) by the 
detached vertex cuts. The background composition is illustrated in Fig. [TT]a. The 
branching fraction measured by ALEPH, B(b -> 57) = (3.11±0.80±0.72) x 10~ 4 , 
is consistent with the CLEO measurement and the Standard Model predictions. 

The other LEP experiments were not able to detect a b — ► 57 signal and set 
upper limits consistent with the CLEO and ALEPH measurements: 
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Figure 10: Inclusive EL spectra in the CLEO b — ► 57 measurement obtained with 
the event-shape analysis (left) and with inclusive B reconstruction (right), (a) 
T(4S) data (solid histogram), scaled below-thre°hold data (dashed histogram) 
plus estimated T(45) backgrounds (points with en or bars), (b) Background- 
subtracted data (points) and Monte Carlo prediction based on Ref. |5 
shape of the b — ► 57 signal (solid curve). Note that the range of E 
different for the left and right plots. 
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Figure 11: Inclusive E* spectrum in the ALEPH b — > S7 measurement. On the 
left: data (points), total estimated background (solid histogram), ir° background 
from b — > c decays (dotted histogram), final state radiation background (shaded 
histogram), and all other backgrounds (dashed histogram). The latter comes 
mainly from ir° decays from non-6 sources and from rj decays. On the right: 



background-subtracted data (points) and Monte Carlo prediction [51| for the 
shape of the 6^57 signal (solid histogram) . 
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Figure 12: The X s mass distribution in B — > X s "f decays from the CLEO inclusive 
B reconstruction analysis. The solid curve is a fit to the expected distribution 
from a spectator model. The dashed curve shows the non-/T*(892) component of 
the fit. 



DELPHI II < 5.4 x 10~ 4 , L3 < 12 x 1(T 4 (90% C.L.) 
Combining the CLEO and the ALEPH results, we obtain: 

B(b si) = (2.54 ± 0.57) x 10~ 4 



4.2.2 Theoretical implications 

Dividing the measured value of B(b —* sj) by the Standard Model predictions, 
we obtain: 



V* 



tb 



0.85 ± O.lO(experiment) ± 0.04(theory) , 



consistent with the unitarity constraint [13]: 



V* 

Vt 



Vcb 



lb 



1.01 ±0.18 



Using the measured values to eliminate V t b = 0.99 ± 0.15 (5f| and V c ^ = 0.040 ± 
0.002 H|, we obtain: 

\V U \ = 0.034 ±0.007 

The agreement between the measured and the Standard Model rates (including 
the CKM matrix unitarity) leaves little room for non-standard contributions. 
For example, limits on anomalous WW^ couplings can be obtained [^7| . Such 
anomalous couplings can arise from internal structure of the gauge bosons or 
loop corrections involving new particles. They are parameterized by Ak and A, 
which are both zero in the Standard Model. Non-zero values of these parameters 
would generate anomalous magnetic dipole and electric quadrupole moments of 
the W: m = ^(2 + Ak + A), Q e w = -j^O- + Ak ~ A )- The y would 
either increase or decrease the b — > 57 rate. The region of the Ak — A space 
consistent with the CLEO measurement is shown in Fig. O. The region allowed 
by a pp — > X measurement by the DO experiment pqj is also shown. The 
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two types of measurements are complementary. The b — > s 7 measurement alone 
excludes the U(l) theory in which neutral bosons of electromagnetic and weak 
interactions do not mix. 




-3-2-10 1 2 



Figure 13: Limits on the anomalous WWj coupling parameters A and Ak. The 
hatched regions are consistent at 95% confidence level with the b — ► s 7 rate mea- 
sured by CLEO. The leading-log calculations, together with their uncertainties 
were used to obtain these bands J5(| . The yield of pp — » Wj X measured by DO 
limits the parameters to the interior of the ellipse (95% C.L.) f5"8fl . The Standard 
Model (broken U(l)xSU(2)) and pure U(l) theory are displayed. 



A charged Higgs boson can be exchanged instead of the charged W in the 
penguin loop. Models with two Higgs doublets are divided into two categories: 
Model I - both up-type and down-type quarks get their masses from the same 
Higgs doublet; Model II - up- type quarks get masses from one Higgs doublet, 
whereas down-type quarks get masses from the other Higgs doublet. The free 
parameters of these models are the mass of the exchanged Higgs and tan /3, which 
is the ratio of vacuum expectation values for the two Higgs doublets. In Model 
I the 6^57 rate is decreased relative to the Standard Model prediction. Since 
the CLEO measurement is somewhat below the Standard Model expectation, the 
data are consistent with Model I and a small Higgs mass. In Model II, the Higgs 
contribution always adds constructively to the Standard Model rate. Recently, 
next-to-leading order QCD corrections have been calculated for this model [O]. 
Using both the CLEO and ALEPH measurements, we obtain B(b — > s^)/B(b — > 
civ) < 3.4 x 10~ 3 (95% C.L.) which translates to a lower limit on the charged 
Higgs mass with slight tan/3 dependence: Mh > 490 GeV for tan/3 = 2, as 
illustrated in Fig. 14. The limit is almost tan/3 independent for tan/3 > 2. 

Minimal supersymmetric extensions of the Standard Model (MSSM) include 
Model II charged Higgs doublets. However, the above lower limit on Mh does 
not directly apply to the supersymmetric model since a chargino-stop loop may 
destructively interfere with the charged Higgs and W-top contributions [^] . Nev- 
ertheless, the data impose interesting constraints on the minimal supersymmetry. 
Naively speaking, either charged Higgs, chargino and stop are all heavy or all 
light. 

The limits on new physics imposed by the measurements can be expressed in a 
model independent way by bounds imposed on Wilson coefficients. The 6 — > s 7 
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Figure 14: Exclusion region in the plane B(b — > sj)/B(b — * civ) versus Mjj for 
tan/3 = 1 (dashed line), tan/? = 2 (solid line), and tan/3 = 5 (dot-dashed line) as 
calculated by M. Ciuchini et al. [jl4| including next-to-leading QCD corrections 
and combining theoretical uncertainties linearly. The upper limit on B(b — > 
s^f)/B(b — > civ) is indicated (horizontal line), together with a corresponding 
lower limit on Mh for tan /3 = 2 (vertical line) . 



decay rate is sensitive mostly to the value of the Ci coefficient, with slight C§ de- 



pendence as illustrated in Fig. [15| |6(J. This figure also illustrates possible MSSM 
models, with additional theoretical symmetries motivated by supergravity. All 
of these models are made consistent with all direct searches for supersymmetric 
particles at LEP and the Tevatron. They are also out of reach of LEP-II. We see 
that the b — > 57 results severely constrain the minimal supergravity models. 

For constraints on other extensions of the Standard Model imposed by the 
measured 6^57 rate see e.g. Ref. [61]. 



4-3 b — > g?7 and V t d 



2. 



Detection of b — > d 7 is difficult because the rates are suppressed by \Vtd\ /\Vts, 
B(B -> dj) = (0.017 - 0.074) x B{B -» a 7) ||. Rejection of the dominant 
background from b — > S7 decays requires good particle identification, except 
for the simplest exclusive final states in which kinematic cuts alone are very 
effective. CLEO searched for B — > (p, u) 7 decays [47|. No evidence for the signal 
was found due to lack of sufficient experimental statistics (2.6 x 10 6 BB pairs). 
The following upper limits were set (90% C.L.): B(B° -> p° 7) < 3.9 x 10~ 5 , 
B(B° ->w 7 )<1.3x 10~ 5 , and B(B~ -> p~ 7) < 1.1 x 10~ 5 . The ratio B(B -> 
(p, cj) j)/B(B — > 7) can be used to determine | V^| 2 /| Vt s | 2 after corrections 
for phase space and SU(3)— flavor symmetry-breaking effects. Unfortunately the 
latter are somewhat model dependent. Long distance interactions may further 
complicate the analysis |6^, |j~3[ |. From the present experimental limits CLEO 
obtains: \V td \ 2 /\V ts \ 2 < 0.45-0.56, where the range indicates the uncertainty in 
the theoretical factors. 
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Figure 15: Bounds on the contributions from new physics to C-j and C%. The 
region allowed by the CLEO data corresponds to the area inside the two bands. 
Various minimal supergravity models consistent with the direct searches for su- 
persymmetric particles at LEP and the Tevatron are represented by points pop. 



4.4 b^sl+l- 

In the Standard Model, the b — > sl + l~ decay rate is expected to be nearly two 
orders of magnitude lower than the rate for 6^57 decays pC|] . Nevertheless, 
the b — ► sl + l~ process has received considerable attention since it offers a deeper 
insight into the effective Hamiltonian describing FCNC processes in B decays 



[13fl. While b — > 57 is only sensitive to the absolute value of the C7 Wilson 
coefficient in the effective Hamiltonian, b — ► sl + l~ is also sensitive to the sign of 
C7 and to the Cg and C10 coefficients, where the relative contributions vary with 
l + l~ mass. These three coefficients are related to the three different processes 
contributing to b — > sl + l~: b — > 57* — > sl + l~ , & — > ,sZ* — > sl + l~ , and the box 
diagram (see Fig. ||). Processes beyond the Standard Model can alter both the 
magnitude and the sign of the Wilson coefficients. 



4.4.1 Searches in exclusive modes 

The simplest allowed final states are B — » Kl + l~, and B — » K*l + l~. Each 
of them is expected to constitute ~ 10% of the total 6 — ► s rate. The 
most sensitive searches for these decays were performed by the CDF and CLEO 
experiments. 

The CDF search [||] is based on 17.8 pb" 1 of data (~ 5 x 10 s bb pairs for 
I77I < 1) and a di-muon trigger. The backgrounds are suppressed by transverse 
momentum cuts (Pt(m) > 2, P t {n 2 ) > 2.5 GeV, P t (K&) > 2 GeV, P t (B) > 6 
GeV), a detached vertex cut (ct(B) > 100/xm), an isolation requirement and 
a B mass cut. The resulting di-muon mass distributions are shown in Fig. pi]. 
Signals from the decays B — ► K^xfjt ) can be seen. Since the branching fractions 
for these decays had been measured previously by other experiments, CDF used 
these signals for normalization. Reconstruction efficiencies are roughly 0.13% 
for the K, and 0.07% for the K* modes. A few events observed outside the 
tp (hatched) and ip' (cross-hatched) bands are consistent with the background 
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estimates. They find the 90% C.L. upper limits B(B~ — > K~ < 1.0 x 10 5 

and B(B° -» K*° (i+/jr) < 2.5 x 10" 5 . 

The CLEO II experiment searched for these decays in a sample of 66 pairs two 
orders of magnitude smaller (~ 2.2 x 10 6 BB) than in the CDF analysis, though 
with efficiencies larger also by two orders of magnitude (~ 15% for K and ~ 5% 
for K*) and suitably low backgrounds. Thus, by coincidence, the sensitivity of the 
CDF and CLEO II experiments were very similar. In addition to the limits in the 
di-muon mode, B(B~ -> K~ fJ, + fx~) < 0.9xl0~ 5 and B(B° -> K*° (J, + fJ,~) < 3.1 x 
10 -5 , CLEO also set limits using di-electrons: B(B — > K e + e ) < 1.2 x 10 5 
and B(B° -» K* e + e~) < 1.6 x 10" 5 , 

The experimental limits are an order of magnitude away from the Standard 
Model predictions. 



4.4.2 Inclusive searches 

The new CLEO analysis [jjJJ looks for inclusive b — > sl + l~ decays using the 
inclusive B reconstruction technique previously described for b — > 57 decays. 



The obtained di-lepton mass spectra are shown in Fig. 17. Again clear signals for 
B — > A,^ and i? — ► X s ^' are observed. Events outside the ip and tp' bands are 
consistent with BB background estimates (the qq background is small). With a 
sample of 3.3 x 10 6 BB pairs and reconstruction efficiencies around 5%, CLEO 
sets 90% C.L. upper limits, B(b -> se + e~) < 5.7 x 10~ 5 and B(b -> s[i + fi~) < 
5.8 x 10~ 5 (combined: B(6 -» sZ+Z") < 4.2 x 10~ 5 ), The SM predictions @, 
B(b -> s e + e~) = (0.8 ± 0.2) x 10" 5 and 5(6 -> s /i + /x~) = (0.6 ± 0.1) x 10" 5 , are 
again an order of magnitude below the experimental limits. 

The upper limit on inclusive b — > s /x + /U~ previously presented by the UAl 
experiment at SppS collider |j(J has been shown recently to be based on over- 
estimated sensitivity ^|. A similar analysis recently completed by the DO 
experiment at the Tevatron resulted in a less stringent limit, < 32 x 10~ 5 [67], 
than achieved by CLEO. 



4-5 b — > svv 

The rate for b — * svv is enhanced compared to the 6 — ► s l + l~ decays primar- 
ily by summing over three neutrino flavors (b — > st + t~ has a small expected 
rate and will be difficult to detect experimentally). The predicted rate is only a 
factor of ten lower than for 6^57 ^T|: (3.8 ± 0.8) x 10~ 5 . In principle, these 
decays are the cleanest theoretically among all penguin decays. Therefore, a mea- 
surement of the inclusive rate for this process would be of considerable interest. 
Unfortunately, the neutrinos escape detection making it difficult for experimen- 
talists to control the backgrounds. So far, only LEP experiments have been able 
to probe these decays by requiring very large missing energy in a hemisphere 
p§| . |69[ ]. Semileptonic backgrounds are reduced by eliminating events with an 
identified lepton in the signal hemisphere. A detached vertex in the opposite 
hemisphere suppresses non-66 backgrounds. The missing energy distribution in 
a 6— hemisphere obtained by ALEPH [68] in a sample of ~ 0.5 x 10 6 66 pairs is 



shown in Fig. [18[ A signal would be seen as an excess of events at large missing 
energy. The lack of such an excess is used by ALEPH to obtain a 90% C.L. limit, 
B{b^svv) < 7.7 x 10- 4 . 

The exclusive mode B — > K* vv is expected to constitute about 30% of the 
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Figure 16: Di-muon mass distributions in the CDF search for b — > s ^ + via 
exclusive final states S + — ► K + [i + [i~ (left) and S° — > K*° [i + n~ (right). 
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total rate 0. DELPHI obtains 90% C.L upper limits §|: B(B° -> K*° vv) < 
1.0 x KT 3 and B(B S -> <f>vv) < 5.4 x 10" 3 . 

The inclusive limit set by ALEPH is an order of magnitude away from the 
expected rate. Unfortunately, no more data are expected at the Z° peak at LEP. 
Perhaps T(4S) experiments will be able to develop analysis techniques which 
can probe these decays with future high statistics data samples. It is hard to 
imagine that experiments at hadronic colliders will ever have any sensitivity to 
these decays. 
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Figure 18: Missing energy in a hemisphere for the selected bb events by ALEPH 
(points). Shaded histograms show the estimated background distribution. The 
expected b — > svv signal shape is indicated by a dotted line. The two highest 
bins are used to set the upper limit. 



4-6 B — > 77 

The vertical penguin B°{B®) — > 77 (Fig. |I]a) decay rate is expected to be of order 

10- 8 (10- 7 ) [ 

rate by ~ 50% W%- The L3 detector at LEP has a high resolution calorimeter, 
which is well suited to this analysis, with a photon resolution better than 2% for 
energies greater than 1 GeV [ff^]. Using their .B-enriched data sample, L3 looks 
for pairs of high momentum photons with invariant mass near the B mass. The 
data are fit simultaneously for B° and B® signals with an exponential background. 
(No candidates are found in a ±2<r mass window.) The L3 90% C.L. upper limits 
are given in Table ||. Unfortunately, these limits are several orders of magnitude 
away from constraining new physics. 




QCD corrections have been found to enhance the B 
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Table 2: Results for vertical penguins. We give the 90% confidence level upper 
limit on the branching fractions (UL B) and the theoretical predictions. 
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for B — > £ + £ are given in Table |. The 



4.7 B^£+£- 

The Standard Model predictions |fr6] 

search for B° — ► e + e~ and B — ► n + n~ at CLEO [73] is similar to their other rare 
B decay searches. The small background is suppressed using mild event shape 
requirements. No signal events are observed in either channel giving limits on 
each channel of 5.9 x IO -6 . 

Since the CDF experiment has excellent muon identification and a muon trig- 
ger, B — > fi + fi~ decays could be recorded with good efficiency [74|. B® mesons 
are also produced at CDF. Background rejection is achieved using a detached 
vertex and an isolation requirement. In the CDF Run IA and IB data sample 
(w 100 pb" 1 ) they find one candidate which falls into both the B° and B2 mass 
windows. This event is consistent with background, but is assumed to be signal 
for calculating upper limits. Results are included in Table ||. 

Searches for B — > l + £~ have also been done by L3 [75]. The small background 
is suppressed by requiring large B energy. They observe no significant signals 
and fit the data simultaneously for B° and B® to extract upper limits on the 
branching fractions. The L3 limits for B — ► e + e~ are 1.4(1.0) x IO -5 , 

less stringent than the CLEO and CDF limits. The L3 limit for B® — ► e + e~ is 
5.4 x 10 -5 , the only limit for this mode. Their limit for B® —> fx + [i~ , 3.8 x 10~ 5 , 
is less restrictive than the CDF limit. 



5 GLUONIC PENGUINS 

Some representative Feynman diagrams relevant for charmless hadronic B decays 



are shown in Fig. [L9[ In addition to the penguin diagrams of primary interest in 
this review, we indicate examples of b — * u tree diagrams in Fig. |Ts|c, d, and g. 
Such diagrams are suppressed for decays with a strange particle in the final state 
but are dominant for many decays with no strange particle. Since penguins play a 
significant role in such decays without strange particles and the two processes can 
be difficult to separate experimentally, we will include results for both processes 
in our discussion. However we do not include results for decays involving the 
dominant b — > c decay mechanism. 

The first indication of signals in charmless hadronic B decays came in a 1993 
CLEO publication j7(| in which a quite significant signal was found for B — > 
h + h~ , where h is either K or ir. Statistics were not sufficient to obtain significant 
signals for B° — ► K + 7r~ or B° — > tt + tt~ separately. Subsequently CLEO updated 
these results [77], still without an observation of either mode individually, and 
provided limits for many related modes. 
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Figure 19: Feynman diagrams for some of the penguin and tree processes which 
are expected to be dominant for the modes described in this paper. 



Several LEP experiments have used the excellent vertex resolution provided 
by their silicon vertex detectors to obtain virtually background-free evidence for 
charmless hadronic B decays. Examples of such events are shown in Fig. || for 
the ALEPH experiment [78]. The LEP experiments also have some ambiguity 
between the decays B° -> K + n~, B° -> vr+vr", and 5° -> K + K~. Fig. || shows 
the mass distribution for ten candidate charmless hadronic B decays from the 
DELPHI experiment [53]. The final states include 7r + 7r~, K + ir~ , K + K~ , /0°7r + , 
K*°tt + , K+p°, and K + a^. 



5. 1 Exclusive decays of B mesons 

In this section we discuss the recent observation of several new decay modes 
by CLEO and many upper limits which are beginning to approach theoretical 
expectations. 



5.1 Exclusive decays of B mesons 
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ALEPH 




Figure 20: Event from the ALEPH experiment showing a B — ► K + tt decay 
candidate as reconstructed with use of the ALEPH silicon vertex detector. 
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Figure 21: Invariant mass distribution obtained by DELPHI for ten charmless 
hadronic B decay candidates, background events at lower mass, and the Monte 
Carlo expectation for the background. 
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5.1.1 Experimental method 

In addition to the 3.3 million BB pairs from T(4S) running (see Table ||), a 
sample of events obtained at a center-of-mass energy below BB threshold is used 
for studies of backgrounds from continuum qq production. Resonance states 
are reconstructed from charged tracks and photons with the decay channels: 
rf — ► 7/7r + 7r~, rf — > p°7, K° via K$ — > 7r + 7r _ , p° — > 7r + 7r~, p + — > 7r + 7r°, 7r° — ► 77, 
»7 -» 77> V vr+vr-TT , w -> vt+tt-tt , K+if-, K*° -> K+vr-, #*+ -> K+vr°, 
and -► J PC°vr + . 

The primary means of identification of B meson candidates is through their 
measured mass and energy. The dominant background process for all decays 
considered here is continuum e + e~ — * qq production. Signal events have a total 
energy consistent with the beam energy (5290 MeV), while most background 
events have smaller energy. As discussed in Sec. |4.l| , CLEO uses the variable 
AE = -Ecand — ^bearn, which peaks at or near zero for signal events. Signal 
events also have mass M consistent with mg (5280 MeV), only 10 MeV below 
the maximum value, while most background events have much smaller mass. For 
states decaying to the resonances mentioned above, the resonance mass is also 
important in discriminating against background processes. In the case of vector- 
pseudoscalar decays and the rf — > p°"f channel, the helicity angle distribution is 
also used. For modes with a high-momentum charged track or a charged track 
paired with a ir°, dE/dx information is also used to identify the charged track as 
a pion or kaon. 

The large qq background can be reduced by about an order of magnitude with 
the use of event shape information. Since B mesons are produced nearly at rest, 
BB events tend to be spherical, whereas qq events tend to be quite collimated 
("jetty"). 

In order to keep the efficiency high while still effectively rejecting qq back- 
ground, a maximum likelihood (ML) fit has been performed for all of the recent 
CLEO analyses. The inputs to the fit are the quantities discussed in the previ- 
ous two paragraphs, while the outputs are the number of signal and background 
events. This procedure provides an efficiency of 20-50% for most modes, at least 
a factor of two larger than previous CLEO and ARGUS analyses, with a compa- 
rable effective background. 



5.1.2 B — > Kir and related decays 

The simplest modes, both experimentally and theoretically are the decays to two- 
body final states without resonances. This includes Kir and irir final states with 
both charged and neutral pions. The CLEO results |7^] include observation of 
the decay B° — ► K + ir~ , strong evidence for the decay B + — » K°ir + , and limits 
for the other five final states. The signals are summarized in Table |3| the limits 
for the penguin modes are summarized in Table § while those for the modes 
dominated by b — > u processes are included in Table ||. 

Fig. ^ shows the likelihood contours for the three cases with significance 
greater than three standard deviations. While the significance for the K + ir° 
and 7r + 7r° final states are both below 3a, there is strong evidence for their sum 



(hir°). This is similar to the case of the h + h~ final state several years ago [76] 
Fig. |23| shows projections of the fit onto the M and AE axes; for all projection 
plots, cuts have been made on other ML variables to better reflect the background 



5.1 Exclusive decays of B mesons 



29 



near the signal region. 



Table 3: Results for the five decay modes which have signals with significance 
larger than 3a. The second efficiency column includes B r , the product branching 
fraction of resonance secondary decays. 



Signal Theory 

Decay mode events Signif. e(%) eB r {%) g ( 10 ' 5 ) ff(10~ 5 ) 

B u -> K+tt- 21.6l£g 5.6o-t 44 44 1.5+^ ±0.1 0.1-2.4 

B + -> K°ir + 9.2l|j 3.2o-t 35 12 2.3±^;J ± 0.4 0.5-2.0 

B+^ V 'K + 33.1±|;| 7.5a 28 15 6.5l};| ± 0.9 0.7-4.1 

B°^r]'K° 7.ll|: J 3.8a 25 4 4.7±§J ± 0.9 0.9-3.3 

B + ^uK + 12.2±|J 3.9o- 28 25 1.5l§;% ± 0.2 0.1-0.7 

f The significance does not include systematic uncertainties. 

Note that all of the decay modes summarized in Table ^ some of which will 
be discussed in the next section, are decays which are expected to be dominated 
by penguins. There is still no direct evidence for the b —* u tree modes, though 
signals with significance of more than two standard deviations in the ir + ir~ and 
7T + 7r° channels suggest that such decays may be observed soon. This dominance 
of penguins over b — > u channels suggests that the problem of penguin pollution in 
CP violation measurements for modes such as B° —* ir + ir~ may be particularly 
severe. 

While few conclusions can be reached with the present level of precision in the 



Kir system, Fleischer and Dighe, Gronau and Rosner [33] have pointed out the 
promise of such decays for measurements of the weak phases a and 7. This is of 
great importance since these are the two phases that will be hardest to measure 
with time-dependent asymmetry techniques at BaBar and Belle. Subsequently 
Fleischer and Mannel f92| pointed out that if the ratio 



B(B°^K+ir-)+B(B°^K-ir+) 

B(B+ K°ir+) + B(B- KH~) [ ' 

is significantly smaller than 1, useful bounds can be placed on the weak phase 7, 
namely sin 2 7 < R. From Table ||, we find R = 0.65 ± 0.38, indeed less than 1 
but certainly not significantly so. Improvements upon the original concept have 
recently been suggested p3| , |94]]. There are potential complications in this analy- 
sis due to SU(3) symmetry-breaking effects, electroweak penguins and especially 
final-state interactions as discussed recently by many authors @, ||. The sit- 



uation is not yet resolved but it appears promising that future measurements of 
these decay rates, with samples 10-100 times larger than existing samples, may 
be able to constrain the CP phase 7. 

There is considerable interest in the CP asymmetry for these and other decay 
modes. However, since the observed samples are so small, experimentalists cannot 
yet make meaningful measurements. The asymmetry, A = , where N + 

and N_ are the number of K + and K~ events respectively, is a measure of direct 
CP violation. For most modes, predictions are A ~ 0.1, while experimental 
errors are presently greater than 0.3. 
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Figure 22: Likelihood contours for: (a) B° — > K + ir and B° — > 7r + 7r ; (b) 
-> ^+vr° and 5+ vr+vr ; (c) 5+ K iT+ and 5+ -> £T°7r+. 
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Table 4: Results from measurements of penguin-candidate B decay modes. We 
give the 90% confidence level upper limit on the branching fractions (UL B), and 

K* 



the range of theoretical predictions [80, |8lL |82|]. K*, K\ and are shorthand 
designations for K*(892), ifi(1400) and K|(1430), respectively. Modes indicated 
by "**" have been observed (see Table |||). Limits below the upper end of the 
theoretical range are indicated in bold face. Experiment key: AR (ARGUS 
CI (CLEO I 
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Table 5: Results from measurements of other charmless B decays. We give the 
90% confidence level upper limit on the branching fractions (UL B), and the 
range of theoretical predictions |8C, 32 1. Limits below the upper end of the 
theoretical range are indicated in bold face. Experiment key: AL (ALEPH |78| ), 
AR (ARGUS pi), CI (CLEO I PJ), C2 (CLEO II |77|, M M, M). 
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5.1.3 B — * r/'K, B — > ojK and related decays 



Decays with resonances tend to be more difficult to observe because more particles 
are involved and the secondary branching fractions to observable final states are 
frequently below 50%. The observation of signals for B — » r(K [|88| and B — » loK 
|89| | was unexpected, since previous predictions were that the branching fractions 
would be smaller than that observed by at least a factor of two. 

CLEO finds a strong signal for B + — ► rj'K + in both the rf — > i]tt + it~ (5.2a) 
and 7/ — ► /3°7 (4.8<r) channels. Combining these with evidence from the chain 
7r + 7r _ 7r° yields a significance (including systematic errors in the 



77 — > rjTT^TT , 7/ 



yield) of 7.5cr as shown in Fig. 24a. The combined significance for the B° — ► t/I^ 
decay is 3.8cj as shown in Fig. 
shown. 



25|a. Projections onto the B mass axis are also 



3300398-007 




#(B + -^VK + )(10 ) 



5.28 

B mass (GeV / c 2 ) 



Figure 24: (a) Likelihood contours and (b) mass projection plot for B + — ► rj'h + . 
For (b), the curves are an overlay of the best fit function (solid) and background 
component (dashed). The histograms show submodes: 7/ — ► 7/7T7r (77 — > 7r + 7r _ 7r°, 
dark shaded), 77' — > ?y7T7r (77 — > 77, light shaded), and 77' — ► 797 (open). 

Similarly, CLEO finds a signal with 3.9a significance for B + — * ujK + as shown 



in Fig. 2£. The mass projection plot is also shown. The results for all three 
signals are summarized in Table ||. As in the previous section, there are hints of 
signals for the corresponding tree decay modes B + — * rj'iT + and B + — > U7T + , but 
the significance of each is only about 2a. 

The recent CLEO work includes limits on many other modes involving 77, 77', u 
and 6 mesons. Tables | and | summarize these and other results. The subset of 
these results from the recent CLEO measurements is also summarized in graphical 



form in Figs. 27 and 28 



A number of features of these decays were predicted in advance. For instance, 
Lipkin pointed out [96] that due to interference between the penguin diagrams 



shown in Fig. |19| e and 19 f, the rate for B — > 77.fr is suppressed and the rate for 
B — > rj'K is enhanced. Detailed predictions by Chau et al. |^] and subsequent 
authors show this pattern, but the predicted enhanced rate for B — > 77'!^ (~1- 
2xl0 -5 ) was still far below the recently measured value. Gronau and Rosner 
suggested [I34J the possibility that there could be significant enhancements to the 
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Figure 25: (a) — 21n£/£ max and (b) mass projection plot for B° — ► rj'K , with 
curves the same as in Fig. |24|. The histograms show sub modes: rj' — ► r/7T7r with 
77 — > 77 (shaded) and 77' — > /ry (open). 
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Figure 26: (a) Likelihood contours and (b) mass projection plot for B + — > u;h + . 
In (b), the curves are the same as for Fig. ^4], with the addition of a dotted line 
showing the B + — > ujK + component of the fit. The histograms show B + — > uiK + 
(shaded) and B + — * ujtt + (open). 
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rate for the decay B — > rj'K due to flavor-singlet (hairpin) amplitudes. 

It is worthwhile commenting on the evolution of the theoretical predictions 
given in Tables || and [|. The basic method of using effective Hamiltonian theory 
has remained unchanged since the original Bauer, Stech and Wirbel [32| calcula- 
tions in 1987. However there have been significant improvements in the detailed 
implementation: effective Wilson coefficients [97]; knowledge of the CKM matrix 
elements p9| ; values of the relevant quark masses and quark mixing angles (es- 
pecially important for r] and rj' decay channels); proper treatment of the QCD 
anomaly |)S| ]; and more complete treatment of non-factorizable contributions us- 
ing the empirical color factor £. Thus the newest predictions 81], epitomized by 
the recent comprehensive calculations of Ali, Kramer and Lii, represent a sub- 
stantial advance over the older predictions (8(J . In order to try to capture a sense 
of this progress, we have indicated the range of newer calculations with solid lines 
in Figs. ^7] and 28, while the range of calculations prior to 1997 is indicated by 
a dashed line. 

Several of the CLEO limits restrict the range of recent predictions, as indicated 
by bold face in the Tables || and |H| and overlap of the "X" and the theory line 
in Figs. ET] and 28. In other cases they are able to eliminate some theoretical 



hypotheses which have been advanced to account for the surprisingly large rate for 
B — > rj'K. For instance the conjecture that the rate is enhanced by a substantial 
cc admixture in the r/', i.e. rj-r\ c mixing flOOj l ; appears unlikely. Such a mechanism 
would also yield a large rate for B — > n'K* , in contradiction with the limit given 
in Table |j. Several authors [101] now agree that the cc admixture in the rf is small 
and actually leads to destructive interference, hence smaller rate, for values of the 
phenomenological parameter £ < 0.3 (see Sec. 2.6), where the rate predictions 
are largest. 

(2-4) xl0~ 5 



26, 



A typical range of predictions for the rate for B — > n'K is 
Some theorists believe that the present experimental value is a fluctuation and 
will end up in this range pq ], while others suggest that the rate can be enhanced 
by a variety of mechanisms: use of an even smaller value of the strange quark mass 
than is now being used [102]; use of somewhat larger values of form factors [ 102 , 



103]; use of smaller values of the CKM phases 7 84] ; including non-factorizable 



contributions arising from the QCD anomaly [104, 105 , 106], or modification to 
the Wilson coefficients whereby the value of the color-factor £ is different for odd 
and even penguin coefficients [83]. In order to account for the B — > n'K data, 
Deshpande et al. |84| used constructive interference for the cc admixture (now 
apparently ruled out) and the CKM phase 7 = 35°; a recent fit @ of all relevant 
data excludes this value: 7 = (64 ± 12)°. Further data should help distinguish 
among the various possibilities. The rate for B — > rj'K* may prove to be crucial 
in this regard. 

The situation for the uj and cj) decays is quite unclear at present. There have 
been several recent theoretical papers specifically addressing vector-pseudoscalar 

|108|| , with an emphasis on the B + — > u>K + , B + — > lott + , 
4>K* decays, where there are measurements or good limits. 
The theoretical estimates for B + — > ojK + |26|, |4[] tend to be <10 -5 except when 
£ ~ or £ > 0.5. The decay B + — > cj)K + is closely related but there is as yet no 
evidence for this decay; the CLEO 90% CL limit of 0.5 x 10~ 5 tends to exclude 
the large £ range. The recent paper of Ali et al. suggests that these decays 
belong to a class where factorization may not work well (hence £ may be different 
than it is for other decays), since the largest Wilson coefficients are not present 
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final states I 
B + -» 4>K+, and£ 
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Branching fraction 



Figure 27: Graphical summary of various recent CLEO measurements of charm- 
less hadronic B decays. Limits are denoted by "X", significant measurements 
by (diamond) points with error bars, and recent (pre-1997) theoretical ranges by 
solid (dashed) lines. 
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Figure 28: Same as the previous figure with modes involving rj, rj' and to mesons. 
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in the amplitudes for these processes. Cheng and Tseng [83] come to similar 
conclusions for quite different reasons. 

There have been a variety of other speculations regarding the rates for these 
charmless processes. Ciuchini et al. have suggested that there could be sub- 
stantial enhancements in certain decays due to "charming penguins", where the 
charm-quark contribution in the loop becomes large due to large non-perturbative 
corrections to the effective-Hamiltonian calculations [107]. They predict that final 



states such as pK, <j)K* , mK, and ojK* would be substantially enhanced, though 
even this model is unable to account for the large measured rate for B + — > uK + . 
Some limits presented here, especially that for B + — > K°p + , tend to indicate 
that such enhancements are not huge. 

5.2 Non-resonant decays 

The processes discussed so far have been two-body or quasi-two-body decays. 
Results also have been reported for a variety of higher multiplicity, non-resonant 
decays. Such studies are important since it is not obvious how much of the final 
state Kirir, for example, is K*tt, Kp, non-resonant, or some other possibilities. 
Additional incentive for measuring such decays was provided by a paper []109([ 
pointing out that interference effects in the 3tt final state could lead to CP viola- 
tion, and asymmetry measurements could potentially measure the weak phase 7. 
However, since this topic is unrelated to penguins, we do not discuss it further. 
The current experimental situation for these higher multiplicity decays is sum- 
marized in Table |(| In a few cases the limits are becoming restrictive but most 
are far away from the theoretically expected branching fractions. For instance 



the prediction for B+ -► 2it+tt- is 10-50 xl(T 6 |09]. 



5.3 Decays of B° s mesons 

There have been relatively few studies of the B® meson since the mass, 5370 
MeV, is too large for production at the T(4S) resonance, so CLEO and ARGUS 
have not studied these decays. Thus all results for B® charmless hadronic decays 
are from the various LEP detectors. The physics is very similar for B® as for B® 



mesons; imagine the u or d spectator quark in Fig. |19| replaced by an s quark, 
with appropriate changes to the final state mesons. Thus the analogue of the 
penguin decay B — > Kit is B® — > KK. In Table [5], we summarize the present 
knowledge concerning rare B® decays. The limit for B® — * K + K~ is close to the 
expectation, but all others are not. However, new information on these decays 
can be expected in the future from hadronic colliders. 

5.4 Inclusive decays 

There have been a variety of inclusive searches for gluonic penguins, which typi- 
cally involve final states with (hidden or open) strangeness. Inclusive decays are 
intriguing partly because direct CP violation might be observed as a difference 



between the numbers of positively and negatively charged kaons [115]. We de- 
scribe several searches for inclusive decays, beginning with attempts to measure 
the inclusive b — > sg* rate and including one quite surprising observation. 

Even though the gluonic penguin b — > sg* does not have a good signature, 
several experiments have used ingenuity to try to measure the inclusive process. 
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Table 6: Results from measurements of non-resonant B decay modes. We give 
the 90% confidence level upper limit on the branching fractions (UL B). Modes 
above the line are expected to be dominated by penguins. Experiment key: AL 
(ALEPH pOU), AR (ARGUS pip , C2 (CLEO II D (DELPHI g|). 



B° final 
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B+ final 
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state 
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AR 
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AL 
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AR 
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AL 
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AL 
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AL 
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AR 
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AR 


11000 









Table 7: Results from measurements of B® decay modes. We give the 90% 
confidence level upper limit on the branching fractions (UL B), and the range of 
theoretical predictions [27, 113 1 . 



Bg decay mode Experiment UL B (10 6 ) Theory (10 6 ) 
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K+K- 
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pp 

7T 7r° 
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DELPHI 
ALEPH 
ALEPH 
ALEPH |78] 
L3 




46 
210 
170 
59 
210 
1000 
1500 



3-21 
10-18 



0.001-0.02 
0.4-5.6 
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The ARGUS collaboration searched for b —* sg* in samples where one B of a BB 
event is fully or partially reconstructed. They searched for decays of the other B 
involving a kaon and multiple pions, a typical signature of the b — * sg* process. 
They found two events where no possible sub-combination is consistent with any 
charmed particle and set a 90% CL upper limit of 8% on the branching fraction 
for b — ► sg* |11G ], 

DELPHI |117| looked for b — ► sg* by searching for an excess of high pt kaons, 



since more energy is available for kaons from b — > s decays than from b — > c — > s 
decays. A fit to the kaon pt distribution provided a preliminary limit of B(b — > 
sg*) < 5% at 95% confidence level. 

There have been two indirect searches for b — > sg* . CLEO set an upper limit of 



6.8% [ |118|| by accounting for all other types of B decays. DELPHI [119] performs 
a similar search, and finds a limit of 4% at 90% confidence level [120]. While none 
of these searches is yet sufficiently sensitive to observe the ~ 1% signal expected 



from Standard Model processes, they appear to exclude models |38|, |40j, |4l|, |4 
where the rate is enhanced by an order of magnitude (see Sec. 2. 



CLEO has searched |12l| [ for a <j> meson accompanied by an X$ system con- 
sisting of a charged or neutral K meson and zero to four pions, of which at most 
one can be a ir°. The X$ system was required to have a mass less than 2 GeV, 
corresponding to a <fi meson momentum of ~2.1 GeV/c; in this region the back- 
ground is small for <j> mesons arising from b — > c processes, potentially allowing 
detection of a b —* sg* signal. CLEO observed no signal and, using the model of 
Deshpande et al. (l22|, they set an upper limit B(B <pX s ) < 1.3 x 10~ 4 . The 
theoretical expectation for this process is (0.6 - 2.0) x 10" 4 jl2H. 



The only positive evidence for inclusive charmless hadronic B decays is from 
the CLEO analysis of the decay B — > rj'Xs [123]. The technique is the same as 
for the <() inclusive search except that the kaons were required to be charged. 
The momentum of rj' mesons, reconstructed with the decay chain rj' — > r]ii + ir~ , 
rj — ► 77, was required to be in the range 2.0 < p v i < 2.7 GeV/c in order to reduce 
background from b — > c processes. The values of AE and M, as defined above, 
were required to satisfy \AE\ < 0.1 GeV and M > 5.275 GeV. 



The rj mass distribution is shown in Fig. 29; a clear signal of 39 ± 10 events 
is seen for on-resonance data and none for the below-threshold sample. The 
signal, obtained by subtracting the scaled off-resonance data in bins of X$ mass, 
is plotted in Fig. |30|. Note the four events corresponding to B + — ► r]'K + and the 
absence of events in the i<C*(892) mass region, both consistent with the exclusive 



results discussed in Sec. 5.1.3j . Also shown in Fig. 3C are distributions for 



potential background modes such as B — > Drj and B — > D*rj'. Though these 
also tend to have large X$ mass, they are more peaked than the data. These 
and other studies suggest that the observed signal does not arise primarily from 
color-suppressed b — > c decays, though it is difficult to rule this out completely 
without better models of such processes. The efficiency was calculated assuming 
that the signal arises solely from gluonic penguin decays, with an equal admixture 
of Xs states from the kaon up to K|(2200). The efficiency of (5.5 ± 0.3)% leads 
to B{B -> i]'X s ) = (6.2 ± 1.6 ± 1.3) x 10~ 4 for 2.0 < p n > < 2.7 GeV/c. The 
systematic error is dominated by the uncertainty in the Xs modeling. 

Many theorists suggest that the B — > rj'Xs result is the most surprising of those 
included in this review — the theoretical expectation for the high-momentum rj' 
branching fraction is ~ 1 x 10~ 4 . Atwood and Soni [104] first suggested that 
the very large rate could be due to the QCD anomaly via b — > sg* — > srj g. 
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Figure 29: The rj' mass distribution for (a) on-resonance and (b) below-threshold 
data. 
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Figure 30: Distribution of Xs mass for data (points with error bars) and possible 
backgrounds: B — > Drj (solid histogram) and B —>■ D*r{ (dashed histogram). 
The normalization of the backgrounds is arbitrary. 
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Later Fritzsch [124] suggested a similar anomaly-mediated b — > sg process. Both 
processes have a hard X$ mass spectrum in rough agreement with that shown 
in Fig. |30|, but Fritzsch's 2-body process would be falling by 2.5 GeV, while the 
3-body decay of Atwood and Soni would still be rising. There have been several 
other recent papers discussing B — > rfXs 1 102 , 102, 125, 126]. They consider the 
anomaly and the possibility of non-SM contributions. It is still not clear whether 
one of these anomaly models can quantitatively account for the B — ► rj'Xg result, 
though there don't seem to be any better explanations at present without invoking 
new physics. Since other inclusive processes such as B — * <pX$ are not enhanced, 
it seems likely that this result is indeed an "anomaly" involving the rf meson. 

In a similar search for high-momentum r\ mesons, CLEO finds no evidence for 
a signal and sets a 90% CL upper limit for the branching fraction of 4.4 x 10~ 4 . 
The theoretical expectation is that the rate for r] mesons would be suppressed 
relative by about an order of magnitude since their flavor-singlet component is 
small. 
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Table 8: Comparison of future B experiments. Parameters which do not change 
between different experiments at the same collider are entered only once. 
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0.2 
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PP 
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500000 


75000 
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5 ■ 


10" 3 



Atlas 500 
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1 Asymmetric beam energies. * Forward detector. 



Future measurements require more produced B mesons and good background 
rejection. The CLEO II experiment has a substantial amount of new data in the 
analysis pipeline (see Table |]) and continues to accumulate even more statistics. 
There are a variety of new experiments poised to do rare-S physics in the near 
and longer term as summarized in Table p|. 

The CESR machine will be upgraded to higher luminosity and the CLEO de- 
tector will get a particle identification device and a new tracking system (CLEO 
III detector [127]). Also under construction are two new T(45) colliders, the 
Positron-Electron Project-II (PEP-II) at the Stanford Linear Accelerator Cen- 
ter (SLAC) with the BaBar |128|] experiment and KEK-B at the high energy 



laboratory in Japan with the Belle experiment [129]. These new machines have 
asymmetric beam energies for indirect CP violation measurements, which will 
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also help background suppression because of the detached B vertex (see the 
/?7ct column in Table |8|). All "B Factory" experiments (BaBar, Belle, CLEO 
III) are expected to come online in 1999. They should be able to reach branching 
fractions of 0(1O -7 ) at the design luminosity. Unlike CESR, the new colliders 
have a double ring structure, thus in principle, they offer better potential for 
higher luminosity. There are some ideas of how to convert CESR into a double 
ring symmetric collider (CESR IV) for luminosity upgrades beyond the CLEO 
III phase fl30f . 

The clean experimental environment at e + e~ colliders is ideal for background 
suppression when looking for penguin processes which have very small branching 
fractions. Unfortunately, the e + e~ — > bb cross-section is only 1 nb, making it 
very difficult to reach extremely rare decay modes. Hadronic colliders, with bb 
cross-sections approaching 1 mb, may offer the ultimate experimental avenue to 
the secrets of B physics. On the other hand, the total cross-section at hadronic 
colliders is much larger, making background rejection the dominant experimental 
obstacle. The detector recording rate is limited by technological constraints. B 
data can easily be lost by the inability of an experiment to reduce background 
rates to manageable recording rates. It should be emphasized that hadronic 
colliders are built with the primary goal of studying ultra-high-energy processes 
which produce particles with large transverse momenta to the beam {Ptj- In 
contrast, B decays produce particles with transverse momenta not much higher 
than ordinary beam interactions. So far, high-P^ detectors at the Tevatron (CDF 
and DO) have been able to trigger on B physics only in di-muon modes. Since 
off-line background suppression has been accomplished by detached B vertex 
cuts, inclusion of a detached vertex requirement in the trigger may be the best 
triggering strategy for non-leptonic rare decays. 

There have been some dedicated experiments to study c and b quark physics 
with hadronic beams in fixed-target mode in which a hadronic beam (tt~ or p) col- 
lides with a stationary nuclear target (N). While the fixed target charm program 
has played a complementary role to e + e~ —* cc studies, the fixed target experi- 
ments barely succeeded in observing some inclusive signals from tree-level b —* c 
decays of b quarks | 131 |. Since most of the incident beam energy is wasted pro- 
viding the motion of the center-of-mass, the effective collision energy was rather 
low resulting in a very small bb cross-section. To compound these experimental 
difficulties, all beam and target fragments entered the detection apparatus to- 
gether with B decay products. The latest attempt to study B physics in fixed 
target mode is represented by the HERA-B experiment [132] which will collide 
the proton beam from the Hadron-Electron Ring Accelerator facility (HERA) at 
DESY with a wire target inserted into the beam pipe. Use of modern technologies 
offers hope for success, though as in the previous fixed target experiments, the 
HERA-B experiment will have to contend with cross sections as low as in e + e~ 
colliders and with hadronic backgrounds as high as in pp colliders. The first run 
is foreseen for 1998, with a full capacity run in 1999. 

The Tevatron collider is in the midst of an upgrade program with turn-on 
planned for 2000. Luminosity will be substantially increased, with further lu- 
minosity upgrades possible in the long term future ("TeV33"). Improved vertex 
detectors in CDF-II Jl33|, [134| and the DO upgrade |13§ will help background 
suppression. The DO experiment will acquire a magnetic field, but the tracking 
system will remain very compact. CDF-II will implement a B — > h + h~ trigger. 
Unfortunately, even if this trigger is successful, CDF-II will have difficulty distin- 
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guishing B° —* 7r + 7r~ from B° — ► iT + 7r~ and -Bl? — > K + K~ because of the lack 
of high momentum particle identification. 

In the farther future (2005), the LEP tunnel will house a new pp collider, the 
Large Hadron Collider (LHC), with high P t experiments ATLAS |136|] and CMS 
[137]. The larger center-of-mass energy will increase the bb cross-section and 
improve the signal-to-background ratio. Since the event rate at the LHC design 
luminosity will be too high to record bb data, the B physics program in these 
high P t experiments may be limited to initial lower luminosity running. 

Both the Tevatron and LHC programs contemplate installation of dedicated b 



experiments (BTeV [|138| 1 and LHC-B | 139 |). The machine luminosity would be 
reduced for these experiments to keep data rates at manageable levels. Without 
the constraints of high Pj physics, the detectors can be optimized for b physics. 
Both detectors would operate in the forward rather than the central region. For- 
ward detector geometry offers enough space for efficient particle identification in 
the entire momentum range (see below). Furthermore, the decay length of B 
mesons is longer than it is in the central region, allowing for better background 
suppression in the off-line analysis and in the trigger. Last but not least, the 
entire trigger bandwidth can be saturated with b physics events. 

Improvement in particle identification is also crucial at e + e~ machines. The 
separation of gluonic b —* s penguin decays from tree-level b — > u decays is 
already marg inal for CLEO II (Sec. |T|) and there are more ambitious goals. 
We would like to have measurements of b — > d penguin decays despite large 
backgrounds from the dominant b —* s decay; measurement of the electroweak 
penguin B + — » (f)ir + with a large B° — ► (pK + background will be a tremendous 
challenge. The BaBar experiment will use a Detector of Internally Reflected 
Cerenkov (DIRC) to measure the Cerenkov angle of light from particles traversing 
quartz bars. The DIRC gives at least 2.5a tt-K separation up to the maximum 
momentum of 4 GeV. The Belle experiment uses an aerogel threshold Cerenkov 
detector to identify kaons with an 8% misidentification rate for momenta between 
1 and 3 GeV. The CLEO III detector will use a Ring Imaging Cerenkov (RICH) 
detector with a solid radiator to achieve 4<r tt-K separation at a momentum of 2.6 
GeV/c, the maximum value from B decays. HERA-B will have a gaseous RICH 
detector, which is expected to provide 90% efficient kaon identification with less 
than 2% misidentification background. CDF-II and DO are limited to kinematic 
separation of kaons and pions, with little additional help from dE/dx. The BTeV 
and LHC-B designs include gaseous RICH detectors. 
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The penguin program already has been a profitable one, but much penguin 
physics is still left to be done. The measured b — > s'j rate has been found 
to be consistent with Standard Model predictions and has provided interesting 
constraints on new physics models. Since these measurements are still statistics 
limited, more data will provide improved constraints and even allow measurement 
of \ Vtd\ from b — > penguin decays. 

Limits on the rates of electroweak penguins b — ► s£ + £~ and vertical penguins 
such as b — > £ + £~ are currently orders of magnitude above the predictions of 
the Standard Model and give no evidence for non-SM effects. These processes 
need to be measured to confirm Standard Model predictions or to distinguish 
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between different types of new physics. Some penguin modes such as b — > svv or 
B° —* 77 will be very difficult to observe, since they are not suitable for detection 
at hadronic colliders and the e + e~ machines are likely to be limited by the bb 
statistics. 

Several suspected gluonic penguins have been observed. The B — > Kir rates 
are consistent with theoretical predictions. However, the B — ► u>K, B — ► j/if 
and 5 — ► r/'X s rates are all surprisingly large! These rates have stimulated 
many Standard Model and non-SM ideas and more data are needed to resolve 
the situation. More work needs to be done to separate the gluonic penguin 
amplitudes from electroweak penguin or b — > u amplitudes in hadronic decays. 
The rates of penguin decays such as B — > Ktt, B° —* (j)K^ and B° — ► rfKg must 
be measured precisely to give us insight into CP violation. The measurement 
of many penguin modes is needed to give us a handle on penguin pollution in 
indirect CP violation. 

The outlook for penguins is very promising. Within the next five years we 
should have good handles on most electromagnetic and gluonic penguin modes, 
and perhaps a beginning of an understanding of suppressed modes which are 
dominated by electroweak penguins or annihilation diagrams. Perhaps there can 
be a reunion in the year 2002 on the 25th birthday of the penguin and we can 
look back at the tremendous progress that has been made! Darts anyone? 
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